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The search for strongly inversion asymmetric topological insulators is an active research field because these 
materials possess distinct properties compared with the inversion symmetric ones. In particular, it is 
desirable to realize a large Rashba spin- splitting (RSS) in such materials, which combined with the 
topological surface states (TSS) could lead to useful spintronics applications. In this report, based on first 
principles calculations, we predict that the heterostructure of BiTeI/Bi 2 Te 3 is a strong topological insulator 
with a giant RSS. The coexistence of TSS and RSS in the current system is native and stable. More 
importantly, we find that both the Z 2 invariants and the Rashba energy can be controlled by engineering the 
layer geometries of the heterostructure, and the Rashba energy can be made even larger than that of bulk 
BiTel. Our work opens a new route for designing topological spintronics devices based on inversion 
asymmetric heterostructures. 



The study of the effects of spin-orbit coupling (SOC) has been one of the central themes in physics research in 
the past decade. A significant discovery is that a strong SOC can convert a normal insulator into a topological 
insulating phase, in which metallic surface states appear in the bulk gap and they are protected by a nontrivial 
bulk topological order 1 ' 2 . Another important consequence of SOC is the Rashba effect in inversion asymmetric 
structures 3 ' 4 . It has been extensively studied in the context of spintronics, which aims to utilize the electron's spin 
degree of freedom for device applications. It is of great interest to search for a material hosting both robust TSS 
and large RSS which could be used as a playground to study the interplay between these two interesting effects and 
also as a base material for technological applications. 

A nonzero RSS requires structural inversion symmetry breaking. It occurs naturally at sample surfaces or 
interfaces. For example, RSS has been observed at metallic surfaces 5 " 7 , ultrathin metal films 8,9 , and semiconductor 
heterostructures 10 . In particular, ARPES experiments on the topological insulator Bi 2 Se 3 have shown that RSS 
exists at surfaces with artificially deposited adsorbates 11 " 14 . Such RSS, however, are typically small and sensitive to 
the surface environment. On the other hand, a large RSS of the bulk bands is possible in polar compounds with 
heavy elements, such as the noncentrosymmetric layered semiconductor BiTel which was recently reported to 
have the largest known RSS strength 15 . However, BiTel is a trivial insulator at ambient condition. It was proposed 
that under pressure BiTel might be turned into a topological insulator 16 . But pressure induced structural trans- 
ition might also occur as suggested by experiments 1718 . 

Motivated by the recent progress in the experimental fabrication of layered heterostructures 19 ' 20 , here we 
predict that the coexistence of robust TSS and giant RSS can be realized in the [ (BiTel) m (Bi 2 Te 3 ) n ] heterostructure 
in its native state without any need of strain or pressure (m and n represent the number of layers of each 
constituent in a unit cell). Moreover, we find that both the topological order and the RSS energy can be controlled 
by varying the number of layers. In particular, the Z 2 topological invariants of heterostructure with even (odd) m 
are given by [ 1 ;000] ( [ 1 ;00 1 ] ) . And the RSS strength is tunable by varying the ratio between m and n. Especially the 
structures [(BiTeI) 2 (Bi 2 Te3) 1 ] and [(BiTeI) 3 (Bi 2 Te 3 )!] have Rashba energies even larger than that of BiTel. The 
TSS at different crystalline surfaces are analyzed in detail, which demonstrate characteristic features of a strongly 
inversion asymmetric topological insulator. Our findings thus identify a promising topological material for future 
spintronics and quantum computation applications. 
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Results 

Basic crystal structure. The well known topological insulator material 
Bi 2 Te 3 has a layered structure with five atomic layers (Te-Bi-Te-Bi-Te) 
as a quintuple layer (QL) unit 21 . BiTel has a similar layered structure 
with three atomic layers (Te-Bi-I) as a triple layer (TL) unit 15 . The 
interactions between QLs (or TLs) are of weak van-der-Waals (vdW) 
type, while the interactions within one QL (or TL) are of strong 
covalent or ionic type. Since in-plane lattice constants of Bi 2 Te 3 
(4.383 A) and BiTel (4.339 A) are very close, their heterostructure 
[(BiTeI) m (Bi 2 Te 3 )J, which consists of alternating m TLs of BiTel 
and n QLs of Bi 2 Te 3 , is possible to be epitaxially grown along the 
crystallographic c-axis 22 . The simplest structure [(BiTel) i(Bi 2 Te 3 )i] 
(hereafter referred to as Bi 3 Te 4 I) is shown in Fig. 1(a). In the 
following, we shall choose it as a representative example for the discus- 
sion of general structural and electronic properties of this family of 
heterostructures. 

The optimized lattice parameters and bond lengths without and 
with vdW corrections are listed in Table 1. As can be seen from the 
table, the inter-layer bond lengths (I-Te! and Te 3 -Te 4 ) are much 
larger than the intra-layer bond lengths in both cases. With vdW 
correction, the intra-layer bond lengths are hardly changed, while the 
inter-layer bond lengths are significantly reduced, confirming that 
the interactions between TL and its adjacent QLs are of weak and 
non-local vdW type. The Te 3 -Te 4 and I-Tei bond lengths relaxed by 
vdW-DF method are similar to the inter-layer Te-Te bond length in 
Bi 2 Te 3 (3.748 A) and the I-Te bond length in BiTel (3.924 A) 
respectively, which is reasonable. The dynamic stability of the struc- 
ture is further investigated through the phonon spectrum calcula- 
tions. The calculated phonon spectrum with SOC is shown in 
Fig. 1(c). The absence of imaginary frequencies demonstrates that 
our structure of Bi 3 Te 4 I is dynamically stable. In the following elec- 
tronic calculations, the optimized structure with vdW corrections is 
adopted. 

Giant Rashba spin splitting. The calculated band structure of 
Bi 3 Te 4 I without SOC is shown in Fig. 2(a). As we can see, Bi 3 Te 4 I 
is a direct band gap semiconductor with ~ 0.52 eV, much smaller 
than that of BiTel (—1.2 eV) 23 . The direct band gap is located at A- 
point of the Brillouin zone (labeled in Fig. 1(b)). The states near the 



(a) 
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Figure 1 | (a) Crystal structure of the Bi 3 Te 4 I superlattice consisting of one 
QL of Bi 2 Te 3 and one TL of BiTel along the crystallographic c-axis. There 
are eight nonequivalent atoms in one unit cell with the stacking order of 
Bi 1 -I-Te 1 -Bi 2 -Te 2 -Bi 3 -Te 3 -Te 4 from bottom to top, as labeled in the figure, 
(b) Hexagonal Brillouin zone with the high-symmetry points. The 
projected surface 2D Brillouin zones are also shown, (c) The calculated 
phonon spectrum of Bi 3 Te 4 I with SOC included. 



Table 1 | Lattice parameters and bond lengths of crystal structures 
optimized using GGA-PBE and vdW-DF methods (see Fig. 1 (a) for 
the atomic labels) 

Bond length (A) 





a (A) 


c(A) 


Bi r l 


I-Te] 


Te r Bi 2 


GGA-PBE 
vdW-DF 


4.431 
4.431 


18.555 
16.921 


3.291 
3.295 


4.595 
3.850 


3.082 
3.083 




Bi2-Te2 


Te 2 -Bi 3 


Bi 3 -Te 3 


Te 3 -Te 4 


Te4-Bii 


GGA-PBE 
vdW-DF 


3.292 
3.280 


3.278 
3.254 


3.087 
3.094 


4.249 
3.746 


3.075 
3.076 



band edges are dominated by Bi-6p z orbitals for the conduction band 
bottom, and by Te-5p z orbitals for the valence band top. 

When SOC is turned on, as shown in Fig. 2(b), the energy bands are 
spin-splitted and at the same time the conduction band (local) minima 
are shifted away from the high symmetry points Y and A, which is 
similar to the case of bulk BiTel. This is the signature of a Rashba type 
spin splitting 15 . The RSS is characterized by the Rashba energy E R , the 
momentum offset k 0 and the Rashba coupling parameter oc R = 2E R /k 0 . 
These parameters are illustrated in Fig. 2(c) for the conduction band 
bottom near T point. The RSS at T point is slightly anisotropic. Along 
T-K direction, we have E R ~ 76 meV and ot R ~ 2.54 eV'A, while 
along the T-M direction, the corresponding values are smaller, with E R 
~ 74 meV and a R ~ 2.48 eV* A. The RSS of conduction band near A 
point has a slightly smaller Rashba energy (E R ~ 56 meV), but the 
value of a R (—2.51 eV* A) is close to that of the T point. The value of 
RSS in the heterostructure Bi 3 Te 4 I is much larger than the RSS of 
conventional semiconductor heterostructures (on the order of 
meV). This giant RSS can be viewed as being inherited from the 
BiTel component. The bulk BiTel has a RSS of 116 meV (k 0 ~ 
0.05 A" 1 ) in our calculation, which is in good agreement with earlier 
calculation (113 meV) 23 and experimental ARPES data (100 meV) 15 . 

Topological insulating phase. In addition to RSS, the band order is 
also inverted after the introduction of SOC. In the presence of SOC, 
spin and orbital angular momentum are mixed and the p z states are 
transformed into j = 1/2 states. As a consequence, the conduction 
band bottom, which is mainly occupied by the Bi-p z orbital, is pushed 
downward. Meanwhile the valence band top, which is occupied by 
the Te-p z orbital, is pushed upward. Due to the large atomic SOC 
from Bi and Te, the energy gap around A point is inverted. To 
facilitate our understanding of this band inversion process, we 
gradually increase the strength of SOC (A so ) from zero to its full 
value, and monitor the evolution of the band structure. The local 
band gap at k z = 0.5 plane (around A points) and at k z = 0.0 plane 
(around T point) is plotted in Fig. 2(d). As k S o increases, the local gap 
around Y points decreases from 0.85 eV to 0.12 eV without closing. 
In contrast, the local gap around A points decreases, closes, reopens, 
and finally reaches a value of 0.22 eV. This band inversion induced 
by SOC signals a topological phase transition in the current system. 

To pin down the topological phase of the heterostructure more 
precisely, we calculate the topological Z 2 invariants. Since the inver- 
sion symmetry is explicitly broken in Bi 3 Te 4 I, the parity criterion can 
not apply 24 . Instead, we evaluate the Z 2 invariants v 0 ; (viV 2 v 3 ) in 
terms of the Berry gauge potential and Berry curvature by using 
our recently developed general method within the full-potential lin- 
earized augmented plane wave framework 25 . The calculated 
Z 2 = 1; (001) for Bi 3 Te 4 I is consistent with the previous band inver- 
sion analysis. At first sight, this topological insulating phase seems to 
be inherited from the Bi 2 Te 3 component. However, it should be 
noted that the Z 2 invariants of Bi 3 Te 4 I is different from that of 
Bi 2 Te 3 (Z 2 = 1; (000)). Instead it is more like bulk BiTel under pres- 
sure 16 , which has the same topological invariants. The difference in 
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Figure 2 | The band structure of Bi 3 Te 4 I without (a) and with (b) SOC. For (a), the Bii~p z (Bi in TLs) and Bi 2 ,3— (Bi in QLs) orbital projected 
character of bands in the dashed box is shown on the right side. The size of circle denotes the weight of projection, (c) Band dispersion in the vicinity of 
r points showing the Rashba splitting. The parameters of Rashba splitting are also illustrated in the figure, (d) The local band gap at k z = 0.5 plane 
(around A points) and k z = 0.0 plane (around T point) as a function of the strength of SOC, where X S o ~ 1 represents the intrinsic value of SOC. 
(e) Surface energy spectrum of side-cut surface (ac plane) of Bi 3 Te 4 I. 



the index v 3 reflects the different band inversion point in the 
Brillouin zone. For Bi 2 Te 3 , band inversion occurs at T point, while 
for Bi 3 Te 4 I band inversion occurs at A point. The nonzero value of v 3 
can have interesting physical consequences. For example, it has been 
shown that a system with v 3 = 1 may offer a topologically protected 
one- dimensional conducting channel along the c-axis 26 . 

The nontrivial Z 2 value is tied to the existence of TSS 27 . As another 
verification, we calculate the surface energy spectrum directly. 
Fig. 2(e) shows the local density of states (LDOS) of a side-cut surface 
(the ac plane of Bi 3 Te 4 I). A Dirac cone in the bulk energy gap can be 
clearly seen, which is located at A point of the projected two-dimen- 
sional Brillouin zone. This again confirms the topological insulating 
phase of Bi 3 Te 4 I. 

Controllability of properties. Thus far, we have shown that the 
simplest heterostructure Bi 3 Te 4 I was a strong topological insulator 
with giant bulk RSS. Now we extend our investigation to the family of 
heterostructures [(BiTeI) m (Bi 2 Te 3 ) n ] with different number of layers 
m and n, and see how the topological invariants and the RSS depend 
on the layer configuration. The calculated results are summarized in 
Table 2. We observe that the coexistence of strong topological 
insulating phase and giant bulk RSS is common to all these 
heterostructures. Quite interestingly, the v 3 index of Z 2 has a 
dependence on the parity of m. Namely v 3 = 1 when m is odd and 
v 3 = 0 when m is even. For odd m, band inversion occurs at A point, 
while for even m, it occurs at T point. The location of the surface 
Dirac cone is associated with the band inversion point, hence can also 



Table 2 | Energy gap (E g in meV), Rashba splitting energy [Er in 
meV) and Z 2 invariants in heterostructures with different layer 
configurations. (TL m QL n refers to the structure [(BiTel) m (Bi 2 Te3) n ]) 

TLiQL] TLtQI.2 TI^QLa TL 2 QLi TL 2 QL 2 TL3QL1 

E g 117 1 14 84 22 52 14 

E R 76 38 22 130 62 139 

Z 2 1;(001) 1;(001) 1;(001) 1 ;(000) 1 ;(000) 1;(001) 



be controlled by the parity of m. (Note that this is meaningful for side 
surfaces where the projections of A and T are distinguishable.) From 
the results, we see that the Rashba energy is also tunable to a great 
extent by varying m and n. E R generally increases with larger m and 
smaller n. It should be noted that heterostructures [(BiTeI) 2 (Bi 2 Te 3 )i] 
and [(BiTeI) 3 (Bi 2 Te 3 )!] have extremely large E R of 130 meV (ot R ~ 
5.23 eV'A) and 139 meV (oc R ~ 5.45 eV*A) respectively, which are 
both larger than that of bulk BiTel. 

Surface states. A salient feature of an inversion asymmetric topolo- 
gical insulator is that different crystal surfaces have different TSS. In 
our heterostructure system, due to the giant RSS and the freedom to 
control the properties by layer configurations, more interesting 
physics are expected to appear on the sample surface. As two repre- 
sentative examples with different Z 2 invariants, we shall discuss 
Bi 3 Te 4 I with Z 2 = l; (001) and Bi 4 Te 5 I 2 ([(BiTeI) 2 (Bi 2 Te 3 )i]) with 
Z 2 = 1; (000) in the following. 

Due to the lack of inversion symmetry, there are two different 
surfaces of Bi 3 Te 4 I (or Bi 4 Te 5 I 2 ) perpendicular to the c-axis, namely 
the I -terminated surface and the Te- terminated surface. They are 
constructed by cleaving the adjacent QL and TL (TL and TL for 
Bi 4 Te 5 I 2 ) through the weakly coupled Te-I gap. As can be seen in 
Fig. 3(a,c), although with different v 3 index in Z 2 , the two hetero- 
structures have similar surface spectra. There are two different kinds 
of electronic states in the energy window around the bulk band gap. 
The states with dispersions identical for both I -terminated and Te- 
terminated surfaces are the bulk RSS states whose wavefunctions 
spread over the bulk region, thus are insensitive to the surface envir- 
onment. The other states with linear dispersions that cross the bulk 
gap are the TSS. Owing to the strong inversion asymmetry, TSS are 
very different at different atom-terminated surfaces. For both het- 
erostructures, the Dirac point of TSS stays in bulk band gap at I- 
terminated surface [Fig. 3(a, c)], while it is buried in bulk valence 
bands at Te-terminated surface [Fig. 3(b, d)] . If Fermi level lies in the 
bulk gap below the Dirac point of I -terminated surface, the two 
surfaces of a sample would have different types of charge carriers, 
which can be used to construct a topological p-n junction 28 . 
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Figure 3 | The surface energy spectrums of I-terminated (a), Te-terminated (b) surfaces of Bi 3 Te 4 I; and I-terminated (c), Te-terminated (d) surfaces of 
Bi 4 Te 5 I 2 . The Fermi level is set to the bulk valence band maximum and dotted lines in (a) and (b) denote various constant energy levels where spin textures 
are imaged. In (e)-(g) the spin-resolved constant energy contours are plotted for the I-terminated surface of Bi 3 Te 4 I, while those for the Te-terminated 
surface are shown in (h)-(j). The red and blue colors denote the positive and negative values of spin projections of S xy S r and S zy respectively. 
The arrows label the spin polarizations of the outermost contour. 



Due to the coexistence of two types of spin-momentum locking — 
TSS and RSS, the surface spin texture is expected to be more com- 
plicated in our system. In Fig. 3(e-g), we show the spin projections of 
S x , S y , and S z for the I-terminated surface of Bi 3 Te 4 I at different 
constant energies. In Fig. 3 (el) at energy 250 meV, there are three 
contours: inner, middle, and outer ones. All of them have the largest 
spin distributions along k x direction while suppressed values along k y 
direction. The outer and inner ones are from the RSS states which 
have opposite spin polarizations, while the middle one is from the 
TSS. When shifting the energy level below the RSS crossing point 
[Fig. 3(e2)], the inner spin contour will change its polarization and 
then switches place with the middle TSS contour and merges with the 
outer one [Fig. 3(e3)], leaving only two spin contours. Below the 
conduction band edge, only the TSS contour remains and its spin 
distribution changes sign when the energy level crosses the Dirac 
point [Fig. 3(e4)]. The spin projection of S y shown in Fig. 3(f) shows 
similar features. Combining results in Fig. 3(e) and 3(f), one can 
easily obtain the in -plane spin chiralities of TSS and RSS states. 
For TSS, the upper and lower Dirac cones have the left- and right- 
handed spin textures respectively. For RSS states, the inner-branch 
has the left- (right-) handed spin textures above (below) the RSS 
crossing point, while the outer-branch always has the right-handed 
spin texture. The S z projection shown in Fig. 3(g) represents the out- 
of-plane spin texture. It can be observed that the outer RSS branch 
has the largest out-of-plane spin component than TSS and the inner 
branch of RSS. For the Te-terminated surface, only the spin projec- 
tions for the upper Dirac cone are shown [Fig. 3(h)-3(j)] because the 
lower Dirac cone is deeply buried in the bulk valence bands. The spin 
textures of RSS states are identical to those of I-terminated surface, 
hence are not shown here. 

Apart from the spin-momentum locking discussed above, recent 
theoretical and experimental works found that in topological insu- 
lator Bi 2 Se 3 there was an additional locking between spin and orbital 
for TSS, leading to a spin-orbital texture 29 " 31 . Following the method 
of Zhang et a/. 30 , we calculate the in-plane (p x and p y ) and out-of- 
plane (p z ) orbital projections together with the orbital -selective spin 



textures for TSS of both I-terminated and Te-terminated surfaces of 
Bi 3 Te 4 I. First let's examine the I-terminated surface. As shown in 
Fig. 4, for the in-plane orbitals, the right-handed spin texture is 
coupled to tangential (radial) orbital texture on the upper (lower) 
Dirac cone, (tangential and radial in-plane orbitals refer to linear 
combinations of p x and p y orbitals with their orbital axis along tan- 
gential or radial directions 30 .) For the out-of-plane p z orbital, it is 
coupled to the left- (right-) handed spin texture on the upper (lower) 
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Figure 4 | The p x and p y orbital- selective spin textures of Bi 3 Te 4 I at the 
upper (a) and lower (b) Dirac cone for I-terminated surface and (c) the 
upper Dirac cone for Te-terminated surface, respectively. The 
corresponding p z orbital- selective spin textures are shown in (d), (e), and 
(f). More red and blue colors in (a)-(c) denote more p x and p y character 
respectively, while more red colors in (d)-(f) denote more p z character. 
The surface states in the energy range of —0.1 eV ~ 0.2 eV are used for the 
orbital projections. The black arrows represent the total in-plane spin 
textures related to the corresponding orbitals. 
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Dirac cone. As for the Te-terminated surface, the orbital -selective 
spin textures are opposite to the I-terminated ones. We have there- 
fore demonstrated that the interesting spin-orbital texture of TSS, 
initially observed in Bi 2 Se 3 3 \ also appears in our inversion asymmet- 
ric heterostructure systems. 

Discussion 

As we mentioned before, in several aspects, the BiTeI/Bi 2 Te 3 hetero- 
structure retains the properties of the bulk BiTel. However, the bulk 
BiTel is a trivial insulator. It can only be converted to a topological 
insulator under pressure. In contrast, for the system which we study, 
no external constraint is required. It is the intercalation by layers of 
Bi 2 Te 3 that changes the topological order of the structure. To better 
understand this result, we start from the orbital analysis of Bi 2 Te 3 
and BiTel. Let's focus on the conduction band edge. The states there 
are mainly from Bip orbitals. For a QL of Bi 2 Te 3 , thep orbitals from 
the two Bi atoms can be combined to form bonding and antibonding 
orbitals. When crystal field is taken into account, p z orbital is split 

from p x>y orbitals, with the bonding orbital p+ = —7= (p B z l2 —p z * 3 ) 

v2 

having the lower energy 21 , as shown in the left shadowed area of 
Fig. 5. Similar energy splitting due to crystal field also occurs for 
BiTel 23 . 

When the layers of two materials are stacked together to form a 
heterostructure, due to the overlap between orbitals, the levels with 
close energies and the same symmetry character are coupled and 
repel each other. Using Bi 3 Te 4 I as an example, on the right hand side 
of Fig. 2(a), we show the Bii~p z (Bi in TLs) andp^ (Bi in QLs) orbital 
projected character of the low energy conduction bands. At T point, 
thep^ and Bii-p z orbital are mainly concentrated in the lowest and 
the second lowest conduction band respectively. As moving towards 
A point, the coupling between p^ and Bi Y -p z orbitals is getting 
stronger, which pushes the lowest conduction band downward. 
Therefore the intercalated Bi 2 Te 3 not only provides the p^~ orbital 
which has a lower energy than the Bii-p z orbital, but also through the 
coupling between the two orbitals it helps to further reduce the band 
gap (from 1.2 eV in bulk BiTel to 0.52 eV in Bi 3 Te 4 I). Hence it 
becomes easier to achieve the band gap inversion when SOC is 
turned on. 

This band gap reduction induced by intercalation also helps to 
explain the observed giant RSS. As been pointed out in the study of 
bulk BiTel 23 , one of the several factors contributing to a large RSS is 
that the band gap should be small. Approximately, the RSS strength is 
inversely proportional to the gap size. The reduced gap size of the 
heterostructure helps to realize a giant RSS. This point is reflected 
from the results in Table 2. The two structures ([(BiTeI) 2 (Bi 2 Te 3 )!] 

QL Bi 2 Te 3 TL BiTel 




Figure 5 | Schematic diagram of the evolution from atomic orbitals of Bi, 
Te and I into the states near the band edges of the BiTeI/Bi 2 Te 3 
heterostructure in the absence of SOC. The evolution stages are explained 
in the main text. 



and [(BiTeI) 3 (Bi 2 Te 3 )!]) with RSS energy greater than that of bulk 
BiTel indeed have very small band gaps. 

To conclude, we have investigated the crystal structure and elec- 
tronic properties of the heterostructures [(BiTeI) m (Bi 2 Te 3 ) J by first 
principles calculations. We discover that these heterostructures host 
both strong topological insulating phase and giant RSS in their native 
states. Moreover, both the Z 2 invariants that characterize the topo- 
logical properties and the Rashba energy that characterizes strength 
of RSS can be controlled by varying the layer numbers m and n. The 
surface electronic structure and spin texture present unique features 
of the coexistence of TSS and RSS, and of strongly inversion asym- 
metric topological insulators 32 . Our result indicates a rich physics in 
this class of heterostructures which provide an ideal platform for 
future development of spintronics and quantum computation 
devices, for example designing the topological p-n junctions, creating 
majorana fermions via combining with superconductors 33 . 

Methods 

Structure optimizations. The lattice parameters and the atomic positions of the 
heterostructures are optimized using the projector augmented wave method 34 
implemented in VASP code 35 ' 36 . Perdew-Burke-Ernzerhof parametrization of the 
generalized gradient approximation (GGA-PBE) is used for the exchange correlation 
potential 37 . Since the interaction between TLs and QLs is of the vdW-type similar to 
bulk BiTel and Bi 2 Te 3 , the vdW corrections by the approach of Dion et a/. 38 (vdW- 
DF) as implemented in VASP is employed to further relax the out-of-plane lattice 
parameter c and atomic positions 39 . The plane wave energy cutoff is set to 300 eV, and 
the Brillouin zone is sampled byal6X 16X4 mesh. Both the crystal's shape and the 
atomic positions are relaxed until force on each ion becomes less than 0.01 eV' A" 1 . 
The phonon spectrum is calculated using the PHONOPY code 40 through the DFPT 
approach 41 . 

Bulk electronic structure calculations. The electronic ground- state calculations are 
performed using full-potential linearized augmented plane- wave (FP-LAPW) 
method 42 , implemented in the package WIEN2k 43 . The convergence criteria are more 
stringent than those of structural optimizations, with i*C max ^MT = 8.0 and a 20 X 20 X 
4 /c-mesh, where £ M t represents the smallest muffin-tin radius (2.5 Bohr for all 
atoms) and K max is the maximum size of reciprocal-lattice vectors. Spin-orbit 
coupling is included by a second-variational procedure 42 . 

Surface electronic structure calculations. Surface electronic structures are calculated 
using the combination of maximally localized wannier function (MLWF) 44 ' 45 and 
surface Green's function methods 46 ' 47 . The construction of MLWF using package 
wannier.90 48 is a non-self-consistent process on a 8 X 8 X 3 /c-mesh with formerly 
converged self-consistent charge potential. The surfaces with different atom- 
terminations are considered as semi-infinite systems, from which the iterative 
method 46 ' 47 is used to obtain the surface Green's function. The LDOS is extracted from 
the imaginary part of the surface Green's function, which is used to analyze surface 
electronic structures. With the advantage of addressing an individual surface, the 
surface Green's function methods is more suitable for the analysis of surface- 
dependent TSS. 
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